use "C:\Users\evanl\Desktop\Projects\COVID\Analysis\Final Cleaned Experimental Data.dta", clear
set seed 271828

gen id=_n

// See R analysis script for excluded observations and implementation of robust standard errors

reg fBehaveW13 test
estimates store a1
reg fBehaveW13 riskaverse
estimates store a2
reg fBehaveW13 i.test##c.riskaverse, hc3
estimates store a3


coefplot (a1, drop(_cons)) ///
			(a2, drop(_cons)) ///
			(a3, drop(_cons riskaverse 1.test)) ///
			, scheme(s1mono) xline(0) ciopts(recast(rcap)) nokey ///
			rename(test = "Treatment" riskaverse = "Risk Aversion" 1.test#c.riskaverse = "Interaction") /// 
			xtitle(Coefficient Magnitude, size(large)) ylabel(, labsize(large)) xlabel(, labsize(medium)) ///
			xscale(range(-0.5 1.5)) xlabel(-0.5(0.25)1.5)

graph export "C:\Users\evanl\Desktop\Projects\Covid\Analysis\Figure 2.png", as(png) replace

			
reg attitudeW13 test
estimates store b1
reg attitudeW13 riskaverse
estimates store b2
reg attitudeW13 i.test##c.riskaverse
estimates store b3


coefplot (b1, drop(_cons)) ///
			(b2, drop(_cons)) ///
			(b3, drop(_cons riskaverse 1.test)) ///
			, scheme(s1mono) xline(0) ciopts(recast(rcap)) nokey ///
			rename(test = "Treatment" riskaverse = "Risk Aversion" 1.test#c.riskaverse = "Interaction") /// 
			xtitle(Coefficient Magnitude, size(large)) ylabel(, labsize(large)) xlabel(, labsize(medium)) ///
			xscale(range(-0.5 1.5)) xlabel(-0.5(0.25)1.5)

graph export "C:\Users\evanl\Desktop\Projects\Covid\Analysis\Figure 3.png", as(png) replace

		
reg perceptionW13 test
estimates store d1
reg perceptionW13 conspir 
estimates store d2
reg perceptionW13 i.test##c.conspir 
estimates store d3


coefplot (d1, drop(_cons)) ///
			(d2, drop(_cons)) ///
			(d3, drop(_cons conspir 1.test)) ///
			, scheme(s1mono) xline(0) ciopts(recast(rcap)) nokey ///
			rename(test = "Treatment" conspir = "Conspiratorial Ideation" 1.test#c.conspir = "Interaction") /// 
			xtitle(Coefficient Magnitude, size(large)) ylabel(, labsize(large)) xlabel(, labsize(medium)) ///
			xscale(range(-0.5 1.5)) xlabel(-0.5(0.25)1.5) 

		
graph export "C:\Users\evanl\Desktop\Projects\Covid\Analysis\Figure 4.png", as(png) replace	
